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Abstract 

We investigate the impact of neutron capture rates near the A = 130 peak on the r-process 
abundance pattern. We show that these capture rates can alter the abundances of individual 
nuclear species, not only in the region of A = 130 peak, but also throughout the abundance 
pattern. We discuss the nonequilibrium processes that produce these abundance changes and 
determine which capture rates have the most significant impact. 

PACS numbers: 26.30.Hj,26.50.+x 



The origin of approximately half of the solar abundances of elements heavier than the Fe 

Soup (Z > 26) is the astrophysical rapid neutron-capture process known as the r-process 
,0). While the basic recipe for generating r-process elements is understood, the production 
of heavy nuclei through repetitions of neutron-capture and /3-decay, much remains to be 
learned about the nuclear physics and the astrophysics of this process. 

On the astrophysics side, a number of environments are under consideration as candidate 
sites. These include core collapse supernova, e.g. 

aaaa 

, neutron star mergers, e.g. 

B 8, black-hole neutron sta r mergers^. *fl, gannna ray bursts, e.g. Q, and shoeked 
surface layers of O-Ne-Mg cores, e.g. [11|, 112] - On the nuclear physics side, there is a lack 
of experimental data for most of the ~ 3000 nuclei that participate in the r-process since 
the majority of these nuclei are highly unstable. With the development of radioactive beam 
technology, an increasing number of these nuclei are now accessible to experiment. 

It is well known that there are two pieces of nuclear data that are important for the r- 

n 

process: beta decay rates and nuclear masses, for a review see [13J. The former has the largest 
leverage on the abundance pattern since the beta decay rates determine the time it takes to 
progress from the seed nuclei to the large A nuclei. Thus in traditional environments, such 
as the neutrino driven wind, the beta decay rates can determine whether an A = 195 peak 
is formed or not. The mass model is generally accepted to be important as well, since the 
equilibrium abundances along isotopic chains are determined by neutron separation energies, 
which are in turn determined by the nuclear masses. Thus, the mass model has considerable 
leverage on the shape of the distribution. 

In contrast, neutron capture rates have received relatively little attention. One might 
be tempted to argue that since the actual values of the cross sections only become impor- 
tant during the last stages of the r-process, i.e. during freeze-out, their impact is limited. 
However, Surman and Engel [ijj] pointed out that neutron-capture rates in the A = 195 
region can influence the r-process abundances in this same region at late times. Addition- 
ally, a number of groups 



la . flf] . [17I. [rsl . Q have examined simultaneous modifications of 

all neutron capture rates; Rauscher [151 showed that such modifications can alter the time 
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until the onset of freeze-out. Beun et al. [20J noted that a single neutron capture rate, that 
of 130 Sn, can affect changes throughout the r-process abundance pattern. In this paper we 
present calculations of the influence of neutron capture rates near the A = 130 region on the 
r-process abundance pattern. We find that the effect is non-negligible, affects nuclei across 



the entire r-process region, and can be comparable to the distribution uncertainties which 
stem from different nuclear mass models. 

Since there is little in the way of measurements for the relevant neutron capture rates, 
uncertainties are hard to quantify. A comparison of three different sets of theoretically 
calculated neutron capture rates [20] shows that most rates vary within a factor of 100, 
although a handful vary up to a factor of 1000 or more. 

The bottom panel of Figure Q] shows two different abundance patterns that we obtain 
by changing the mass model from ETFSI 21] to FRDM [22J]. We use capture rates of 23] 
with the ETFSI masses and of 241 with the FRDM masses; in both cases we use beta decay 

n 



rates from Moller et al. [25|. We use a parameterized neutrino driven wind as in [20j, with 
a timescale r = 0.1 s, a constant entropy per baryon of s/k = 100, and initial electron 
fraction of Y e = 0.26. In the top panel of Fig. [T]we show abundance patterns using the same 
astrophysical conditions, beta decay rates, and the ETFSI mass model, but with the neutron 
capture rates of a range of nuclei with 128 < A < 138 changed by a factor of 100. Since 
abundance plots are usually made on a log scale, differences in nuclear abundances are better 
seen on a linear plot. In the middle panel we show percent differences for both scenarios. 
The average change in abundance is 43% for the neutron capture rates and 52% for the mass 
model, where the latter effect is due primarily to the differing neutron separation energies. 

Neutron capture rates cannot influence the abundance distribution throughout most of 
the r-process when the temperature, T, and neutron number density, Tl n , £1X6 high enough 
for (n,7) ^ (7> n ) equilibrium. This occurs when the neutron capture rate, A n , and the 
photo-dissociation rate, A 7 , between two adjacent nuclei are balanced and the abundance 
distribution along an isotopic chain is determined entirely by the Saha equation. The series 
of most abundant isotopes for each element as determined by the Saha equation is called 
the equilibrium r-process path. We determine the actual abundances and path using a 



dynamical nuclear network calculation 
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201 ] . In the dynamical calculation, each set 



of (n/y) ±=i (7,n) rates falls out of equilibrium at different points in the simulation. A 
comparison between the equilibrium path and the actual path provides a measure of how 
strongly equilibrium is maintained throughout, which we quantify as follows. We calculate 
an abundance averaged neutron separation energy for the system by two different methods: 
a) taking the abundance weighted average of the separation energy, S njactua i, and b) taking 
the equilibrium abundance weighted average of the separation energy, S nieg . When the 



neutron capture and photo-dissociation rates are no longer sufficiently rapid to respond to 
changes in the equilibrium r-process path resulting from changes in the temperature and/or 
neutron density, equilibrium fails, individual rates become important, and S n ^ actua i departs 
from S n)eq , as shown in the top panel of Fig. [2j 

On the left hand side of this panel, the separation energies are rising as the system moves 
toward stability at late times. However the photo-dissociation rates are not rapid enough to 
maintain equilibrium so S niactU ai cannot keep pace with S n%eq . At this point individual photo- 
dissociation rates can become important. In particular, the photo-dissociation of nuclei that 
are highly populated can influence the overall neutron number density. As an example, we 
consider 132 Cd in the baseline simulation of Fig. [2j 132 Cd is along the equilibrium and actual 
r-process path while equilibrium obtains; as S n>eq rises as described above, the equilibrium 
path shifts to 130 Cd. However, the photo-dissociation rate of 132 Cd is too slow to fully shift 
the actual path in the baseline simulation. The photo-dissociation rate of a nucleus (Z, A+l) 
depends on the capture rate of the nucleus (Z, A): 



where (<jv)z,a is the averaged value of the neutron capture cross section, a, and the relative 
velocity, v. We therefore perform a simulation where the capture rate of 131 Cd is increased 
by a factor of 10 and compare with the baseline simulation. With the increased rate, at late 
times 132 Cd is able to photo-dissociate to 131 Cd, which in turn promptly photo-dissociates to 
130 Cd, as depicted in the bottom panel of Figure [2J This late-time shift influences the entire 
abundance distribution as Cadmium is the most abundant element in this simulation — these 
dissociations therefore emit a significant number of neutrons that are subsequently captured 
throughout the nuclear network. In the top panel of Figure [3] we show the number of 
neutrons captured relative to the baseline simulation for both the A = 130 peak and the 
region above the A = 130 peak. The solid line curve decreases sharply, showing that less 
neutrons are captured in the A = 130 peak as compared to the baseline simulation, due to 
the photo-dissociation effect. There is a corresponding increase in the number of neutrons 
captured elsewhere. 

Returning to the top panel of Fig. [2] we see that S n)actua i and S nieg cross at a slightly 
later time. At this point \p > A n , A 7 so the increase in S njactU ai is primarily due to beta 
decay. Nevertheless, some capture and release of neutrons is still possible, with A n > A 7 




(1) 



due to the falling temperature. Therefore, a few nuclei that are highly abundant can drain 
significant numbers of the remaining neutrons from the system. For example, at late times, 
the nucleus 131 Sn becomes highly populated in the baseline simulation because it is on the 
beta decay chain from a waiting point nucleus at late times, as shown in the bottom panel 
of Figure O The neutron capture rate from 131 Sn controls the number of neutrons that are 
captured on the way to 132 Sn. The impact of this drain of neutrons is shown in the bottom 
panel of Figure El This figure shows the number of neutrons captured in a simulation where 
the neutron capture rate of 131 Sn is increased by a factor of 100 relative to the baseline 
simulation. In the A = 130 peak more neutrons are captured and so correspondingly less 
are captured in the region above the A = 130 peak, i.e. the rare earth and A = 195 regions. 



This is the same effect as described in 20 1. 



Changing neutron capture rates creates a situation where more or less neutrons are avail- 
able for capture across the system at different times during the freeze-out process. This 
causes shifts in the final abundance pattern that depend on the relative strength and timing 
of the two effects. We have shown two specific examples, but there are a number of rates 
that can have these same effects. In order to quantify the impact of individual neutron 
capture rates, we calculate an average fractional change, F, by summing the differences of 
all mass fractions, X, of all nuclear species 

F = 100 \XA,baseline ~ (2) 

A 

for simulations in which only one rate is varied. 

Using different nuclear data or different astrophysical conditions can cause shifts in the 
path, i.e. different nuclei are populated during freeze-out from equilibrium. Therefore, 
different neutron capture rates can become important. We study the results from a number 
of different simulations in order to determine the range of rates which can play a role. We 
perform twelve baseline simulations using six different nuclear data sets: 

• FRDM masses and rates, and Moller et al. [3] /3-decay rates 

• FRDM masses and rates, and Moller et al. [25J /3-decay rates 



as above but with experimental data [27| included 



ETFSI masses and rates, and Moller et al. [26] /9-decay rates 



ETFSI masses and rates, and Moller et al. [25] /3-decay rates 



as above but with experimental data 27J] included 



and two different outflow timescales in the wind model, r = 0.1 s and r = 0.3 s. For each of 
the twelve baseline simulations we perform 250 additional simulations in which the neutron 
capture rate of each of the 50 nuclei in Figure H] is varied by a factor of 10, 50, 100, 500, 
or 1000. The darkest color squares in Figure H] represent nuclei for which a factor of 10 in 
the neutron capture cross section yields a change in the overall abundance pattern, F, of 5 
to 20% in at least one of the simulations. Nuclei which require a factor of 50 in order to 
produce at least a 5% change in the abundance pattern have a medium shade. For uncolored 
nuclei, no simulation produced a 5% change in the abundance pattern. Increasing the rates 
of more than one nucleus at a time will compound the effect, producing for example, the 
43% change seen in Figured! Since odd-iV nuclei have smaller separation energies and tend 
to fall out of equilibrium sooner than even-iV nuclei, there is a larger sensitivity to odd-iV 
capture rates as shown in Figure HI 

Studies of astrophysical sites and conditions, future observations, and improved nuclear 
physics input are important for determining the astrophysical site of the r-process. Presently 
the most lacking piece of input is on the side of astrophysics theory. Nevertheless, it is 
advantageous to be able to disentangle which effects arise from astrophysical conditions and 
which arise from nuclear physics considerations. In addition to beta decay rates and nuclear 
masses, neutron capture rates also play a role in determining final r-process abundance 
distributions. 
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FIG. 1: The bottom panel shows two different abundance patterns that were obtained using the 
same astrophysical conditions but different mass models. The top panel shows abundance patterns 
that were obtained using the same mass model, but by increasing fifty neutron capture rates in 
the A = 130 peak region by a factor of 100. The middle panel shows the percent difference in 
abundance for each point on the curves on the bottom panel (stars) and the top panel (squares). 




FIG. 2: Top panel: Shows the equilibrium (dashed line) and actual (solid line) separation energies 
along the r-process path as a function of time for the onset of freezeout in the baseline simulation 
with FRDM masses and capture rates and conditions r = 0.1 s, Y e = 0.26, and s/k = 100. Bottom 
panel: Depicts the nuclear flow in the region of 131 Sn and in the region of 131 Cd for the baseline 
simulation. Shaded squares show the relevant location of the path: at the beginning of freeze-out 
from equilbrium 129 Ag, 132 Cd, and 133 In are highly populated and toward the end 130 Cd and 131 In 
are highly populated. In color version, red arrows indicate photodissociation, blue arrows indicate 
neutron capture, and the black arrow shows beta decay. 
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FIG. 3: Top panel: Shows the effect of increasing the rate of Cd by a factor of 10 relative to 
the baseline simulation. The black line shows the rate of neutrons captured in the A = 130 region 
relative to the baseline simulation and the dotted line shows the rate of neutrons captured above 
the A = 130 region. Bottom panel: Shows the effect of increasing the neutron capture rate of 
131 Sn by a factor of 100 relative to the baseline simulation. More neutrons are captured on 131 Sn 
(and therefore in the A = 130 peak - black line), leaving less neutrons to be captured in the region 
above the A = 130 peak (dotted line). 
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FIG. 4: Shows the nuclei which have neutron capture rates which effect a F ~ 5% or more abun- 
dance change for a rate increase of a factor of 10 (dark color squares), a factor of 50 (medium color 
squares), and a factor of 100 to 1000 (light color squares). 



